On the qubit limit of cavity nonlinear optics 
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Many proposals for solid-state photonic implementations of quantum information processing uti- 
lize high-quality optical resonators to achieve strong coupling between guided fields and hetero- 
geneously incorporated qubits. Given the practical difficulty of accurately placing quantum dots, 
vacancy centers, or other such atom-like emitters throughout a complex nanophotonic circuit, it 
would be natural to consider whether high-quality resonators could be used in conjunction with 
bulk optical nonlinearities to create optically-coupled qubit degrees of freedom via lithographic 
patterning of a homogeneous medium. A recent limit theorem for quantum stochastic differential 
equations can be used to prove rigorously that this should be possible, in principle, using resonators 
incorporating a strongly Kerr-nonlinear, x' 2 '-nonlinear, or two-photon absorbing material with very 
low loss at the fundamental optical wavelength. 

PACS numbers: 03.65.Yz, 42.50.Ex, 42.50.Lc, 42.50.Pq 



Optical resonators incorporating nonlinear refractive or 
absorptive materials have long been of interest in classi- 
cal information technology for photonic computation and 
switching [TJ . Recent advances in the fabrication of high- 
quality nanophotonic resonators via lithographic pattern- 
ing of optically nonlinear materials hold great promise for 
such applications, enabling the study of optical bistabil- 
ity and related switching phenomena at very low (sub- 
femtojoule) energy scales H, H|. In the attojoule (few- 
photon) limit that one might ultimately hope to reach 
with strongly nonlinear low-loss materials and ultra-high 
quality resonator fabrication, cavity nonlinear optics in 
fact crosses over into a regime of nonclassical physics 
characterized by quantum fluctuations and squeezing Q . 
Should we therefore expect some kind of quantum dy- 
namics to take the place of classical (dissipative) switch- 
ing in this limit? In the case of an intracavity Kerr non- 
linearity it can be shown that the two stable branches of 
dispersive optical bistability limit to the two basis states 
of a coherent qubit that responds canonically to optical 
fields driving the resonator. A qubit limit can also in 
principle be reached with nonlinearity and, some- 
what surprisingly, with two-photon absorption as long as 
low losses are maintained at the fundamental wavelength 
of the cavity resonance. It is hoped that this observation 
may help further to motivate research on lithography- 
compatible materials for nonlinear optics. 

Our proof of the above assertions relies on a recent 
limit theorem for quantum stochastic differential equa- 
tions (QSDE's) due to Bouten, Van Handel and Silber- 
farb [5|. We therefore begin by establishing notation re- 
quired to connect with the technical setup of this the- 
orem, however, the reader need not be familiar with 
QSDE's to follow the main results below. The cav- 
ity nonlinear optical models that we wish to consider 
can be specified via choice of coefficients in a Hudson- 
Parthasarathy equation 0, @] for the time evolution of a 
Heisenberg propagator. The coefficients, conventionally 
labeled K, Li and Ny, are generally operator- valued with 
i and j as indices for the n input-output channels cou- 



pled to the cavity. The Hamiltonian of the intracavity 
dynamics is H = Im [K] and the corresponding master 
equation for the intracavity dynamics is given by (H = 1) 
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Note that the Li as defined here and in [j| are hermitian 
conjugates of the usual Lindblad operators from quantum 
optics. Although the Nij do not appear in the master 
equation for the internal degrees of freedom, they are 
required to construct scattering models (Til [l2j in which 
the cavity input-output channels drive or are driven by 
additional quantum systems [j| [T(| ■ 

Within this framework, we first consider a single-mode 
optical resonator incorporating a Kerr nonlinear medium. 
The corresponding Hamiltonian can be written [§[ 

H = A a f a + \ a^a^aa = (A - v + xa t a)o t a, (2) 

where a is the annihilation operator for the intracavity 
field, A is the detuning of the cavity resonance from the 
frequency io p of a rotating frame for the model, and \ is a 
coefficient that fixes the strength of the Kerr nonlinearity. 
The second form of the Hamiltonian emphasizes that a 
Kerr nonlinearity induces a shift in the cavity resonance 
frequency proportional to the intracavity photon num- 
ber. If we consider a cavity that has two input-output 
channels (i.e., is coupled to two waveguides, one of which 
could be used to model distributed losses) with associated 
photon coupling rates K a i,2, we can choose = Sij and 



Li 



(3) 



As can be seen by assembling the corresponding master 
equation, the total cavity decay rate is thus set to n a = 
n a i + n a 2- If we want to consider the dynamics with 
a coherent driving field incident through channel 1, we 
can use the series product [TTl EH to derive the simple 
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FIG. 1: Numerical integrations of the pre-limit (finite \) mas- 
ter equation for a driven Kerr-nonlinear cavity with K a i = 0.5, 
K a2 = A = 0, a = 10 and (a) x = -100, (b) X = -20. Black 
crosses are (0|p|0), red circles are green diamonds are 

the sum of the remaining populations, red solid curve is 
according to the limit (qubit) master equation ©. 

modifications 

H — > An'a + pVaal — (aa' - a'a), 

2i 

Li -> e^VSwat + a*), L 2 e'^V^a 1 - (4) 

Here a is the complex amplitude of the driving field 
( | ck | 2 has units of photons/sec) and we assume that the 
rotating frame of the model has the same frequency 
as the driving field. With this notation in place we 
can remind the reader of the physical origin of disper- 
sive optical bistability in the classical regime, in which 
K a 3> \x\- Assuming A and x have opposite sign and 
that |A| > K a , there can be two equilibrium conditions 
for fixed K a i,2, a, A, x- either (a'a) is small because the 
cavity is driven off-resonance and the detuning 'correc- 
tion' x°^ a is negligible, or (a^a) becomes large enough to 
make (A — \ + X a ^ a ) ^5 K a and the external driving field 
can self-consistently maintain the high intracavity pho- 
ton number because of the decrease in effective detuning. 

If we now contemplate the possibility of achieving 
\x\ 3> Ka in future nanophotonic systems with strongly 
nonlinear media and high-quality resonators, we immedi- 
ately recognize that by setting A = \ m Eq. ([2]) we would 
be left with a model of a highly anharmonic oscillator. If 
the level spacing between the one- and two-photon basis 
states is thus made very different from that of the spacing 
between the zero- and one-photon states, it would seem 
to follow that a finite driving field should drive transitions 
only within the lowest 'qubit' subspace (13(. This type 
of scenario recalls some basic ideas from circuit quantum 
electrodynamics [l4j in which the nonlinear inductance 
of a Josephson Junction can be used to isolate a qubit 
subspace of a superconducting LC-oscillator in a similar 
fashion. 



FIG. 2: Numerical integrations of the pre-limit (finite g, Kb) 
master equation for a driven x^ -nonlinear cavity with K a i — 
0.5, K a2 = A = 0, K b = 5000, a = 10 and (a) g = -3000, 
(b) g = —600. Black crosses are (0a0(,|pj0a0i,), red circles are 
(l a 0&[p|l a 0&), green diamonds are the sum of the remaining 
populations, red solid curve is according to the limit 

(qubit) master equation (HJ). 

To develop this idea more rigorously, we invoke Theo- 
rem 11 from Q to obtain the limit QSDE that governs the 
open quantum system dynamics in the parameter regime 
\x\ ^ K a,a, A. We satisfy the 'singular scaling assump- 
tion' and 'structural requirements' for the theorem by 
setting (in the notation of [f|) 

Y = ixalalaa, A = 0, F 1 =F 2 = 0, 
Ha - span{|0),|l)}, 

OC 

Y = Y , ~\, H(m\, (5) 

where \m) denotes the m-photon Fock state. We thus 
obtain a limit master equation that can be written 

P = ~i[H', p] + K a jo-po-t - ^o-Vp- ^pcrVj, 

H' = Aer^er — i^/K a i(aa' f — ev*er), (6) 

where a = |0)(1| and the evolution of the cavity state is 
now confined to the qubit subspace, span{|0), |1)}. We 
clearly recognize this master equation as being equivalent 
to that of a two-level atom with spontaneous emission 
rate K a , driven by an external field with Rabi frequency 
2-y/Kai \a\ and detuning A. Hence the emergence of qubit- 
like dynamics is established in this limit. Fig. Q] shows 
numerical examples [l5j of the behavior of the pre-limit 
(finite x) master equation with n a \ = 0.5, n a 2 = A = 0, 
a = 10. With x = —100 we observe a damped Rabi 
oscillation within the qubit subspace in close correspon- 
dence with quantitative predictions of the limit model, 
and for \ — — 20 there is still qualitative agreement al- 
though some population leaks quickly out of the qubit 
subspace. 
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A qubit limit can also be derived in the case of an in- 
tracavity x^-type nonlinearity 0| that degenerately cou- 
ples the high-quality cavity mode of interest to a much 
lower-quality pump mode at twice the frequency [l6| . 
Starting from the pre-limit coefficients JVy — 8ij and 



Li 

Li 

H 




(aa f - a* a) + 2A6 f 6 



+i^(aWfe — aab^), 



(7) 



where b is now the annihilation operator for a pump mode 
at twice the frequency of the fundamental, we first con- 
sider the limit of large g and Kb with g 2 /Kb fixed. Note 
that we assume no driving of the pump cavity mode and 
assign it only a single input-output channel for decay. 
Applying once again the theorem from [j| with 
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1 



-Kb 



tfb, A=-(aarf -atatb), 



2 ' 2 
Fi = F 2 = 0, F 3 = e—'-y/Kb 
Ho = span{|0 a & ), |l o 6 ), . . .}, 
-2 
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where I a denotes the identity operator on the Fock space 
of the a mode, we obtain the limit master equation 



-i[H',p} + K 

9 2 J t t 
< aapa'a' 

Kb I 
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pa* a 



—a'a'aap— —pa'a'aa 



H' = Aa^a — i^J K a \{aa) — a* a), 



(9) 



with the b mode now adiabatically eliminated [33j | . We 
can obtain a final master equation identical to Eq. ^ by 
using the theorem again to take the limit of large g 2 / Kb, 
setting 



Y 


5 t t 
= a'a'aa, 


A = 0, 




2Kb 


Fi 


= F 2 = 0, F 3 


= e 2i Up t J_ a t a t 
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= span{0 a o ), 


i.o 6 )}, 
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oo 

= h® V — 

, o rn a 

m a =2 


(m a - 



(10) 



where 7 is the identity operator of the b mode. 

In view of the latter result, we see that we could actu- 
ally start with a model © with (cavity enhanced [17], EH) 
two-photon absorption (TPA) as the only optical non- 
linearity and still expect to obtain qubit-like dynamics 
in the limit of large TPA rate. This is encouraging 
from a practical perspective since nonlinear absorption 
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FIG. 3: (a) Time evolution of {a' a) for the TPA master equa- 
tion @ with K a i = 0.5, K a 2 = 0, A = 0, q = 10 and 
g 2 /K b = 5000 (black), 200 (green), 40 (red), 8 (blue), 0.5 
(black dash-dot). 



seems generally more straightforward to engineer 1 1 91-42 ll 
than nonlinear refraction without concomitant TPA [22j . 
The idea of using strong TPA to confine the quantum 
evolution of the intracavity field state to the zero/one- 
photon subspace is reminiscent of several recent pro- 
posals for quantum logic gates [H, [24[ or all-optical 
switching [l|| [2f| based on the 'quantum Zeno' effect. 
Fig. [5] shows numerical examples generated by integrat- 
ing the pre-limit master equation corresponding to 
the model ([7]); very similar results are obtained by inte- 
grating the pre-limit TPA master equation ([9]) with cor- 
responding values of g 2 /k (not shown) . In both cases the 
apparent decay rate of the Rabi oscillation increases when 
the relevant parameter (g or g 2 /Kb) is not very large, al- 
though the confinement to the qubit subspace remains 
quite good; the manner of approach to the qubit limit 
is thus seen to be somewhat different in the /TPA 
scenarios than for the Kerr nonlinearity. Fig. |3] displays 
plots of (a' a) as a function of time, generated by numeri- 



cal integration of the TPA master equation ([9]) for g 2 / Kb 
decreasing from 5000 down to 0.5, illustrating the loss 
of Rabi oscillations and a transition to monotonic evo- 
lution qualitatively similar to that of a linear cavity. It 
is interesting to note however that with g 2 / Kb — > the 
steady-state photon number with K a = 0.5 and a = 10 
would be (a* a) « 800, hence the optical limiting effect of 
TPA on the response of the fundamental mode to a co- 
herent external drive is seen to be remarkably strong even 
for the lowest TPA rate considered (g 2 / Kb = 0.5 = K a ). 
Indeed, relatively weak (in the context of the current 
discussion) TPA in bulk semiconductors is widely inves- 
tigated in classical photonics for the design of broadband 
optical limiters [26| and related logic devices for pulsed 
signal fields [27| . 

The natural information processing paradigm of purely 
linear photonics is that of Laplace transfer functions act- 
ing on classical analog signals [28j], while the addition of 
moderate optical nonlinearities enables bistable switch- 
ing phenomena that can be used to implement classical 
digital logic [l|. Here we have seen that the incorpora- 
tion of extreme optical nonlinearities (while preserving 
low loss at the fundamental wavelength) could in prin- 
ciple enable a transition to quantum information pro- 
cessing [29j with localized photonic qubits [3(|. It seems 



4 




(b) 
0.15 



(c) 

0.1 




FIG. 4: (a) Time evolution of Sb \p\, a relative-entropy based 
non-Gaussianity measure [3lL [33 |. for the pre-limit TPA 
model with K a i = 0.5, fi; a 2 = 0, A = and g 2 /nb = 200 (red 
solid trace) or g 2 /Kb = 20 (black dash-dot trace), (b) Wigner 
function of the intracavity field at t — 0.2101 for g 2 /Kb = 20. 
(c) Wigner function of the intracavity field for g 2 / Kb = 200. 



natural to ask however if there should be useful 'quan- 
tum analog' signal representations that would naturally 
suit the regime of strong but finite optical nonlinearity, 
and whether a corresponding 'semi-quantum' informa- 



tion processing paradigm could be of practical interest 
as an intermediate stage in the evolution of nanophotonic 
technologies. 

While a serious exploration of the latter idea is beyond 
our current scope, it may be worth noting that strong but 
finite cavity-enhanced TPA could be a useful resource 
for generating non-Gaussian states in nanophotonic cir- 
cuits, which potentially could be used to implement 
ideas from continuous- variable quantum information the- 
ory [3l|. Fig.lUa) displays Sb[p],& relative entropy- based 
measure [3lL [33 ] of the non-Gaussianity of the intracav- 
ity field state as a function of time, computed by nu- 
merical integration of the pre-limit TPA master equa- 
tion © with K a i = 0.5, K a2 = A = and g 2 / n h = 200 
(red solid trace) or <? 2 / K & — 20 (black dash-dot trace), 
starting from |0) at t = 0. The Wigner functions of the 
intracavity field states near the point of maximum non- 
Gaussianity around t s» 0.21 are shown in Figs.[4|b,c). It 
is interesting to note that with g 2 /Kb = 20, for which the 
intracavity state of Fig. [DJb) has nearly 10% population 
in 1 2), the dissipative nonlinearity of TPA can still induce 
a negative component of the Wigner function even in a 
regime of imperfect confinement to the qubit subspace. 
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